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The large N 3sf=4 gauge theory with quenched !N=2 quark matter in the presence of a magnetic field 
displays chiral symmetry breaking. We study the temperature and chemical potential dependence 
of this theory using its gravity dual (based on the D3/D7 brane system). With massless quarks, at 
zero chemical potential, the theory displays a first order thermal transition where chiral symmetry is 
restored and simultaneously the mesons of the theory melt. At zero temperature, these transitions 
with chemical potential are second order and occur at different chemical potential values. Between 
the three there are two tri-critical points, the positions of which we identify. At finite quark mass 
the second order transition for chiral symmetry becomes a cross over and there is a critical point 
at the end of the first order transition, while the meson melting transition remains similar to the 
massless quark case. We track the movement of the critical points as the mass is raised relative to 
the magnetic field. 



D 
> 

in 

oo 

00 

c4 
o 
o 



X 



I. INTRODUCTION 

The phase diagram in the temperature chemical po- 
tential (or density) plane is a matter of great interest in 
both QCD and more widely in gauge theory [1-3]. In 
QCD there is believed to be a transition from a confining 
phase with chiral symmetry breaking at low temperature 
and density to a phase with deconfinement and no chiral 
symmetry breaking at high temperature. In the stan- 
dard theoretical picture for QCD with massless quarks, 
the transition is first order for low temperature but grow- 
ing density, whilst second order at low density and grow- 
ing temperature. The second order transition becomes a 
cross over at finite quark mass. There is a (tri-)critical 
point where the first order transition mutates into the 
(second order) cross over transition. In fact though there 
could still be room in QCD for a more exotic phase dia- 
gram [3] as we will discuss in the context of our results 
in our final section. 

In this paper we will present a precise holographic [4-6] 
determination of the phase diagram in the temperature 
chemical potential plane for a gauge theory that displays 
many of the features of the QCD diagram, although the 
precise details differ. A pictorial comparison of our the- 
ory to QCD can be made by comparing Fig 5 to Fig 10. 

The theory we will consider is the large N N=4 gauge 
theory with quenched 3\T=2 quark matter [7-10] which 
has been widely studied [11]. An immediate difference 
between the D\f=4 glue theory and QCD is that the ther- 
mal phase transition to a deconfmed phase occurs for in- 
finitesimal temperature since the massless theory is con- 
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formal [5]. Essentially the entire temperature chemical 
potential phase diagram of our theory is therefore char- 
acterized by strongly coupled deconfincd glue. 

The quark physics is more subtle though - the phase 
diagram in the temperature chemical potential (density) 
plane for the JsT = 2 quark matter has been studied in [12- 
17]. When the quark mass is zero the theory is conformal 
and the origin of the phase diagram is a special point with 
confined matter. Immediately away from that point, in 
cither temperature or chemical potential, a first order 
transition moves the theory to a deconfined theory (the 
mesons melt [18-20]). 

When a quark mass is present in the 3\T=2 theory the 
meson melting transition occurs away from the origin. 
This transition has been reported as first order with a sec- 
ond order transition point where the first order transition 
line touches the T = chemical potential axis [15, 16] 
(in the grand canonical ensemble). Interestingly there 
is a phase transition line in the temperature versus den- 
sity plane (in the canonical ensemble) in which the quark 
condensate jumps [12, 13]. This area of the phase dia- 
gram is intrinsically unstable though and not realizable 
by imposing any chemical potential [16]. 

The crucial ingredient we will add to the theory is chi- 
ral symmetry breaking which will also bring the theory 
closer in spirit to QCD. As shown in [21-24] the K = 2 
theory in the presence of a magnetic field displays chiral 
symmetry breaking through the generation of a quark 
anti-quark condensate. At zero density the finite tem- 
perature behaviour has been studied [21, 22] and there 
is a first order transition from a chiral symmetry broken 
phase at low temperature to a chiral symmetry restored 
phase at high temperature. In this paper we will include 
chemical potential as well to map out the full phase di- 
agram in the temperature chemical potential plane. We 
will find a chiral symmetry restoration phase transition, 
which is first order for low density and second order for 
low temperature - there is a critical point where these 
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transitions meet. This physics is in addition to a meson 
melting transition which is first order at large temper- 
ature but apparently second order at low temperature. 
This latter region of transition is interesting because it 
is associated with a discontinuous jump from an embed- 
ding off the black hole to one that ends on it and it looks 
naively first order. However, when we plot any avail- 
able order parameter in the boundary theory it appears 
second order. 

We will also track the movement of these transition 
lines and critical points as the quark mass rises relative 
to the magnetic field. The infinite mass limit corresponds 
to the pure ZNT=2 theory without magnetic field [14, 16]. 
The second order chiral symmetry restoration transition 
becomes a cross over the moment a mass is introduced. 
The first order transition structure though remains, even 
to the infinite mass limit, with two critical points: one is 
the end point of the first order transition and the other 
is the the end point of the second order meson melting 
transition. This structure was not reported in the results 
in [14, 16] 1 but this is not surprising since the structure, 
in that limit, is on a very fine scale. We have only found it 
by following the evolution of the larger structure present 
at low quark mass with a magnetic field. In addition we 
present evidence to suggest the parameter space with a 
second order meson melting transition extends away from 
just the T = axis, again, even in the infinite mass limit. 
We have confirmed these results in the strict B = limit 
also. 

The theory we study may appear to be a rather vague 
relative of QCD with magnetic field induced chiral sym- 
metry breaking. On the other hand it is a theory of 
strongly coupled glue with the magnetic field inducing 
conformal symmetry breaking in the same fashion as 
Aqcd in QCD. In fact the magnetic field case in the 
basic N=4 dual is the cleanest known example of chiral 
symmetry breaking in a holographic environment. Other 
deformations of the N=4 gauge theory typically lead to 
an ill-understood IR singular hard wall - see for exam- 
ple [25, 26]. The magnetic field case provides a smooth 
IR wall where we have more control but the results are 
likely to be the same in those more complex cases. We 
can hope to learn some lessons for a wider class of gauge 
theories. 



II. THE HOLOGRAPHIC DESCRIPTION 

The N=4 gauge theory at finite temperature has a 
holographic description in terms of an AdSs black hole 



1 The existence of two critical points is related with the existence 
of the black hole to black hole transition. It is actually just 
visible in Fig 2c of [14] but the authors had not probed it in 
detail previously. After discussion of our results with the authors 
of [14], they have refined their computations and confirmed our 
results. 



geometry (with N D3 brancs at its core) [4-6]. The ge- 
ometry is 



ds 2 



dx 2 



where R 4 = A-Kg s Na 2 and 



/ := 1 - 
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—dr 2 + R 2 dQ 



r H := ttR T 



(1) 
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Here rjj is the position of the black hole horizon which 
is related to the temperature T. 

We will find it useful to make the coordinate transfor- 
mation 



drf_ 

r 2f 



dw 2 



H ' 



(3) 



with wh = th- This change makes the presence of a 
flat 6-plane perpendicular to the horizon manifest. We 
will then write the coordinates in that plane as p and L 
according to 



= vV + L 2 , 



w sm( 



L := wcost 



(4) 



The metric is then 

2 



ds 2 



It'- 



ll 2 

R 2 



(-g t dt + g x dx ) 



— (dp 2 + p 2 dn 2 3 + dL 2 + L 2 dn() , (5) 



where 



(w 4 — wjj) 2 
2w 4 (w 4 + w%) ' 



9t 



H 



2w 4 



(6) 



A. Quarks/D7 brane probes 

Quenched (Nf <C N) N=2 quark supcrficlds can be in- 
cluded in the 3Sf=4 gauge theory through probe D7 branes 
in the geometry [7-10]. The D3-D7 strings are the quarks. 
D7-D7 strings holographically describe mesonic operators 
and their sources. The D7 probe can be described by its 
DBI action 



Sdbi = -T D7 J d 8 ^-dct(P[G] ab + 2ir a 'F ab ) 



(7) 



where P[G] a b is the pullback of the metric and F a t, is 
the gauge field living on the D7 world volume. We will 
use F a b to introduce a constant magnetic field (eg F i2 = 
—F%i = B) [21] and a chemical potential associated with 
baryon number A t (p) ^ [13, 27]. 

We embed the D7 brane in the p and directions 
of the metric but to allow all possible embeddings must 
include a profile L(p) at constant fix. The full DBI action 
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we will consider is then 



S = J d?JCQ>)= (J ea Jdtdxj JdpL(p) , 
where £3 is a volume element on the 3-sphere and 



(8) 



N f T D7 Z- I 1 



W 4 



x^(l + (d p L) 



(w 4 - W%f 



\ 



(9) 



Since the action is independent of A t , there is a conserved 
quantity d (\= -jp-^j and we can use the Lcgendrc trans- 



In the absence of temperature, magnetic field and density 
the regular embeddings arc simply L(p) = m, which is 
the minimum length of a D3-D7 string, allowing us to 
identify it with the quark mass as shown, c should then 
be identified with the quark condensate with the relation 
shown. 

We will classify the D7 brane embeddings by their 
small p behavior. If the D7 brane touches the black 
hole horizon, we call it a black hole embedding, other- 
wise, we call it a Minkowski embedding. We have used 
Mathcmatica to solve the equations of motion for the D7 
embeddings resulting from (13). Typically in what fol- 
lows, we numerically shoot out from the black hole hori- 
zon (for black hole embeddings) or the p = axis (for 
Minkowski embeddings) with Neumann boundary condi- 
tion for a given d. Then by fitting the embedding function 
with (16) at large p we can read off m and c. 



formed action 



S = S- I d^Fpt 



ss 



pt 



63 / dtdx 1 j dp L (p) 



(10) 



B. Thermodynamic potentials 
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£ : 
K 



- N f T D7 {W \..f H) \/K(l + (d P L) I I I ) 



8w 4 
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A2 



(12) 



To simplify the analysis we note that we can use the 
magnetic field value as the intrinsic scale of conformal 
symmetry breaking in the theory - that is we can rescale 
all quantities in (11) by B to give 



L = -N f T D7 (RVB) 4 W ^ H ^ K ( l + L ' 2 ) . ( 13 ) 



/ U! 4 



ur 
1 



(w 4 + w%) 

where the dimensionless variables are defined as 



d 2 ,(14) 



8S 
SA t 



The Hamilton's equations from (10) are d p d 

and d p A t = — |f ■ The first simply means that d is the 
conserved quantity. The second reads as 



doA t 




(17) 



where A* 



\fT2-KU 1 'A t 
RV2B ' 



There is a trivial solution of (17) with d = and con- 
stant At [16]. The embeddings arc then the same as 
those at zero chemical potential. For a finite d, A' t is sin- 
gular at p = and requires a source. In other words the 
electric displacement must end on a charge source. The 
source is the end point of strings stretching between the 
D7 brane and the black hole horizon. The string tension 
pulls the D7 branes to the horizon resulting in black hole 
embeddings [13]. For such an embedding the chemical 
potential(pt) is defined as [13, 27] 



(w,L,p, d) 
I w 



L 



rV2B' i?V2B' RVzb' {Rv r B) 3 NfT D7 2na' 



(15) 



In all cases the embeddings become flat at large p tak- 
ing the form 



~™ ~ c „ 2-Ka'm 
L(p) ~ m + — ,m 



RV2B 



(RV2B) 3 



hm A t (p) 

p—>oo 



dp d 



PH 




(18) 



where we fixed At(pii) = for a well defined At at the 
black hole horizon. 

The Euclidcanizcd on shell bulk action can be inter- 
preted as the thermodynamic potential of the boundary 
field theory. The Grand potential (fl) is associated with 
the action (9) while the Hclmholtz free energy (F) is as- 
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sociated with the Legendre transformed action (10): 

-S 



F(w H ,d) 



N fTDi(R\ -B) 4 Vol 



dp 



"JL^K{l + {L'Y) (19) 



CI(wh,P} := 



-S 



N f T D7 (RVB) 4 Vol 



dp 



-< /(l+(L') a ) 



pa 



K 



W + W 



^ JL ) f? + — P 6 (20) 



where Vol denote the trivial 7-dimcnsional volume in- 
tegral except p space, so the thermodynamic potentials 
defined above are densities, strictly speaking. Since 
K ~ p 6 , both integrals diverge as p 3 at infinity and 
need to be renormalizcd. Thermodynamic potentials, 
(18), (19) and (20) are reduced to B = case if we sim- 
plify omit all -£j and then tildes. See for example (21). 



III. CHIRAL SYMMETRY BREAKING AND 
THE THERMAL PHASE TRANSITION 

We begin by reviewing the results of [21, 22] on mag- 
netic field induced chiral symmetry breaking and the 
thermal phase transition to a phase in which the conden- 
sate vanishes. While they show the embeddings for fixed 
T and different values of B, we will show the embeddings 
for fixed B and different values of T. By fixing B we are 
using it as the intrinsic scale of symmetry breaking in the 
same fashion as Aqcd plays that role in QCD. 

Let us digress here to explain how to understand the 
figures we will present in this paper. For example, in 
Fig 1 we have three columns. The left is the D7 brane 
embedding configuration. The middle shows a plot of 
the allowed values of the condensate c as a function of 
the quark mass fh - these are thermodynamical conjugate 
variables. The right is the corresponding thermodynamic 
potential. Each row is for a fixed parameter we are vary- 
ing - here it's temperature. The left and middle plots 
are plotted by solving the equation of motion (13) with 
the black hole boundary condition that the embedding is 
orthogonal to the horizon. 

The right hand plot is calculated using (19) or (20). 
Both are the same at zero density. We subtract 
limp^oo jp 4 to remove the common infinite component. 

Every point in the middle and right plots corresponds 
to one embedding curve in the left plot. These points are 
color coded with the colors common across each of the 
three plots. The order of colors follows the rainbow from 
the bottom embedding as a mnemonic. 



In the middle plot we can find any transition point by 
a Maxwell construction (an equal area law), which is also 
confirmed by the minimum of the grand potential on the 
right. The vertical dashed line in the middle and right 
hand plots corresponds to the transition point. 

In the left plots the gray region contains embeddings 
that are excluded since they are unstable, as shown in 
the middle and on the right. 

The results for the case of a constant magnetic field 
and varying temperature are displayed in Fig la-d. The 
Fig la (Left) shows the D7 embeddings when T <C B 
and the black hole is small. The embeddings are driven 
away from the origin of the L — p plane - this behaviour 
is a result of the inverse powers of w, when wh <C 1, 
in the Lagrangian (13) which lead the action to grow if 
the D7 approaches the origin (note that the factor of /5 s 
multiplying the action means the action will never actu- 
ally diverge) . There is also a embedding that end on the 
black hole (shown in red) but they are thermodynami- 
cally disfavoured as shown in Fig la (Right). 

At large p the stable embedding with fh — has a 
non-zero derivative so c is non-zero and there is a chi- 
ral condensate i.e. chiral symmetry breaking. The U(l) 
symmetry in the fii direction is clearly broken by any 
particular embedding too. We can numerically read off 
the values of fh and c from the embeddings and their 
values are shown in Fig la (Middle), where the dotted 
blue curves are for Minkowski embeddings, whilst the 
red curves are for black hole embeddings. 

If the temperature is allowed to rise sufficiently then 
the black hole horizon grows to mask the Rl OR. of the 
plane in which the inverse w terms in the Lagrangian are 
large. At a critical value of T the benefit to the fh = 
embedding of curving off the axis becomes disfavoured 
and it instead lies along the p axis - chiral symmetry 
breaking switches off. This first order transition occurs 
at u>h = 0.2516 as shown in Fig lb by Maxwell's con- 
struction (Middle) and by lower grand potential (Right). 
Our value for the critical temeprature agrees with the 

value B = 16 in [22] since our w is the same as in 

[22]. 

We show an example of the embeddings above the crit- 
ical temperature, their grand potential and the evolution 
of the curves in the fh — c plane in Fig lc. 

The Fig Id shows a case when T 3> B when the area of 
the plane in which B is important is totally masked by the 
black hole and the results match those of the usual finite 
T version of the N=2 theory. For fh > uh the embed- 
dings are Minkowski like whilst for small fh they fall into 
the black hole. There is a first order phase transition be- 
tween these two phases which is the meson melting phase 
transition discussed in detail in [28-32]. Minkowski em- 
beddings have a stable mesonic spectrum [10] whilst in 
the case of black hole embeddings the black holes' quasi- 
normal modes induce an imaginary component to the 
meson masses [18, 19]. We can see that the previously 
reported "meson melting" transition at large quark mass 
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(a)Low temperature - wjj = 0.15. Here we see chiral symmetry breaking with the blue embedding thermodynamically preferred 

over the red at rh = 0. 
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(b)Transition temperature - wh = 0.2516. This shows the point where the first order chiral symmetry phase transition occurs 
from the blue to the red embedding. The transition can be identified by considering Maxwell's construction(Middle) or the 

lowest free energy (Right). 
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(c)Above the transition - wjj = 0.3. This is the chiral restored phase with the m = curve lying along the p axis (red). 
w„ = 300 
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(d)High temperature - wh = 300. Here the magnetic field is negligible and the embeddings show the usual finite temperature 

meson melting transition. 



FIG. 1: The D7 brane embeddings(Left), their corresponding m — c diagrams(Middle), and the Free energies (Right) in the 
presence of a magnetic field at finite temperature. (Parameters are scaled or B — 1/2R 2 in terms of parameters without tilde.) 



becomes also the chiral symmetry restoring transition at 
zero quark mass. 



IV. FINITE DENSITY OR CHEMICAL 
POTENTIAL AT ZERO TEMPERATURE 

We can now turn to the inclusion of finite density or 
chemical potential in the theory with magnetic field. In 
this section we consider the zero temperature (uh = 0) 
theory only, and will continue to finite temperature in 
the next section. 



A finite density (chemical potential) at zero tempera- 
ture has been studied in the Jsf=2 theory without a mag- 
netic field in [15], where analytic solutions for both a 
black hole like embedding and a Minkowski embedding 
have been found. When a magnetic field is turned on, an- 
alytic solutions arc not available any more, but we have 
found numerical solutions that continuously deform from 
the known analytic solutions at zero magnetic field. 

Minkowski embedding solutions correspond to zero 
density and finite chemical potential. The black hole like 
embedding is the embedding deformed by the density - a 
spike forms from the D7 down to the origin of the L — p 



(a)Low density - d = 0.01. Here we see chiral symmetry breaking (the blue embedding is preferred over the red embedding) and 

a spiral structure in the fh vs c plane. 
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(b)Increasing density below the transition - d = 0.1. There is still chiral symmetry breaking here with the orange embedding 
preferred to the red. Note the spiral structure in the fh — c plane has disappeared. 






(c) Transition point 
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0.3197. This shows the point where the second order chiral symmetry phase transition occurs. 
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(d)High density d = 1. This is the chiral restored phase with the m = curve lying along the p axis. For larger m the usual 

spike like embedding can be seen. 

FIG. 2: The D7 brane embeddings(Left), their corresponding fh — c diagrams(Middle), and the Free energies (Right) in the 
presence of a magnetic field at finite density. (Parameters are scaled or B = 1/2R 2 in terms of parameters without tilde.) 



plane (Fig 2d (Left)) which has been interpreted as an 
even distribution of strings (i.e. quarks) forming in the 
vacuum of the gauge theory. 

First of all it will be interesting to see how the repul- 
sion from the origin induced by a magnetic field and the 
attraction to the origin by the density compete. Thus 
we start with the canonical ensemble (that is solutions 
with non-zero d) and consider black hole like embeddings 
exclusively. The plot in Fig 2a (Left) shows the embed- 



dings for a small value of density. The solutions show 
the chiral symmetry breaking behaviour induced by the 
magnetic field but then spike to the origin by the den- 
sity at small p. For fh = one should compare the blue 
and red embeddings - the blue one is thermodynamically 
preferred as shown in Fig 2a (Right) The theory shows 
similar behavior to that seen at zero density: there is a 
spiral structure in the fh vs c plane (Fig 2a (Middle)) [21]. 
That will disappear as the density increases. 
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R' H = 




(a)Zero temperature 





wh = 0. The second order meson melting transition and then the second order chiral restoration 
transition are apparent. 



w H = 0.15 






(b)Low temperature - wh = 0.15. The zero temperature structure remains. 






(c)Above the first tri-critical point - wh = 0.23. The meson melting transitions remains second order but the chiral symmetry 

restoration transition is first order. 




(d)Above the second tri-critical point 





wh = 0.24. There is now only a single first order transition for meson melting and 
chiral symmetry restoration. 



w H = 0.2516 
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(e)High temperature - wh = 0.2516. The ground state preserves chiral symmetry for all values of ft. 



FIG. 3: The D7 brane embeddings (Left), their corresponding d — Jl diagrams (Middle), and the grand potentials (Right) for 
massless quarks in the presence of a magnetic field at a variety of temperatures that represent slices through the phase diagram 
Fig 5. (Parameters are scaled or B = 1/2R 2 in terms of parameters without tilde.) 



8 



As the density increases the value of the condensate for 
the fh — embeddings falls - we show a sequence of plots 
for growing d in Fig 2b-d (Middle). There is a critical 
value of d = 0.3197 where c becomes zero for the massless 
embeddings - above this value of d, D7 embedding is flat 
and lies along the p axis (2c-d (Left)) . One can see from 
the plots that there is a second order phase transition to a 
phase with no chiral condensate. In Fig 2c (Left) and 2d 
(Left) we show embeddings at the critical value of d and 
above it respectively. At very large density the solutions 
become the usual spike embeddings of the N = 2 theory 
at zero magnetic field. 

We are not yet done though since there are also 
Minkowski embedding with zero density but constant 
chemical potential. These can have lower energy and be 
the preferred vacuum at a given value of chemical poten- 
tial - that is, they are important in the Grand Canonical 
Ensemble. The relevant analysis is in Fig 3a (Fig 3b- 
c will be explained in the next section). On the left it 
shows the three possible types of embedding of the D7 
for a given chemical potential at zero temperature. The 
black curve is the Minkowski embedding (with d = 0), 
the blue the chiral symmetry breaking spike embedding 
(with d 0) and the red the chiral symmetry preserving 
black hole embedding (with d ^ 0). Strictly speaking 
there is a fourth embedding which lies along the p axis 
and has constant A t = p - its energy is equal for all p to 
that of the red embedding at p = and is never preferred 
over the red embedding with density, so we will ignore it 
hence forth. The trajectory of the three key embeddings 
in the d — p space is shown in the middle plot (note 
that again these two variables are thermodynamical con- 
jugate variables). Finally on the right the grand poten- 
tial is computed. Clearly at low chemical potential the 
Minkowski embedding is preferred and d = 0. There is 
a critical value of p — 0.470 at which a transition occurs 
to the spike embedding. This transition looks naively 
first order since it is a transition between a Minkowski 
embedding and a black hole embedding. However, we 
can see that the Grand Potential appears smooth and 
the quark density is continuous, which is shown again in 
Fig. 4a. The solid lines in Fig 4a are calculated from (18), 
which is based on the holographic dictionary. The dot- 
ted lines are obtained by numerically differentiating the 
grand potential (d = — ^2 ) , which comes from a thermo- 
dynamic relation. This is a nontrivial consistency check 
of the holographic thermodynamics as well as our calcu- 
lation [27, 33]. 

Further in Fig 4b we plot the behaviour of the quark 
condensate through this transition. The density and 
quark condensate are both smooth and the transition 
looks clearly second order. Here we have tested the 
smoothness numerically at better than the 1% level. 
Whether there is some other order parameter that dis- 
plays a discontinuity is unclear but it would be surprising 
that any order parameter were smooth, were the transi- 
tion first order. We conclude the transition is second 




(a) Density: the solid lines are calculated 
from (18) and the dotted lines are 
obtained by numerically differentiating 

the grand potential (d — — §£)• 




(b) The quark condensate. 

FIG. 4: Plots of the order parameters vs chemical potential 
at zero temperature and finite B. Both are continuous across 
the Minkowski to spiky embedding transition(£t ~ 0.47). The 
green arrows indicate the changes of phase. 



order (or so weakly first order that it can be treated as 
second order). This second order nature of the transi- 
tion from a Minkowski to a spiky embedding has been 
shown also in the B = 0,m^0 case at zero temperature 
analytically [15] and numerically [16]. 

Finally, above the chemical potential corresponding to 
the meson melting transition (p = 0.470), non-zero den- 
sity is present and the physics already described in the 
Canonical Ensemble occurs, which turns out to be equiv- 
alent to the results from the current Grand Canonical 
Ensemble. Both Ensemble predict the second order tran- 
sition to the flat embedding at the same point, p = 0.708 
or d = 0.3197, which is the chiral symmetry restoration 
point. Notice that for the Canonical Ensemble we used 
(m,c) conjugate variables on constant d slices, while for 
the Grand Canonical Ensemble we used (p,d) conjugate 
variables on constant fh = slices. This agreement from 
different approaches is another consistency check of our 
calculation. 

On the gauge theory side of the dual, the description 
is as follows. At zero density there is a theory with chiral 
symmetry breaking and bound mesons. As the chemi- 
cal potential is increased d remains zero and the quark 
condensate remains unchanged. Then there is a second 
order transition to finite density (to a spike like embed- 
ding) which is presumably associated with meson melt- 
ing induced by the medium. At a higher density there 



is then a further second order transition to a phase with 
zero quark condensate. 

Finally we note a recent paper [34] that proposed an 
alternative ground state for a chiral symmetry breaking 
theory at finite density. They proposed that the string 
spike might end on a wrapped D5 brane baryon vertex in 
the centre of the geometry. We have not considered that 
possibility here but it might be interesting to investigate 
this in the future. The magnetic field induced chiral sym- 
metry breaking provides a system in which this could be 
cleanly computed without the worries of the hard wall 
present in that geometry. 



0.1 0.2 0.3 0.4 




(a) Wh = 0.15 



0.05 1 0.10 




(b) w H = 0.23 



FIG. 6: Quark condensate vs chemical potential at finite B. 
Both are continuous across the Minkowski (black) to black 
hole (orange) embedding transition. At wh = 0.23 the black 
hole (orange) to black hole(red) transition is discontinous. 



V. THE PHASE DIAGRAM IN THE GRAND 
CANONICAL ENSEMBLE 

We have identified a first order phase transition from a 
chiral symmetry breaking phase with meson bound states 
to a chirally symmetric phase with melted mesons in our 
massless theory in the presence of a magnetic field with 
increasing pure temperature. On the finite density axis 
the meson melting transition is second order and sepa- 
rate from another second order chiral symmetry restoring 
phase transition. Clearly there must be at least one crit- 
ical point in the temperature chemical potential phase 
diagram. We display the phase diagram of the massless 
theory, which we will discuss the computation of, in Fig 
5. 

To construct the phase diagram we have plotted slices 
at fixed temperature and varying chemical potential. We 
display the results in Fig 3a-e where we show the em- 
beddings (Left) relevant at different temperatures, their 
trajectories in the d — Jl plane (Middle) and the grand 
potential (Right). 

The phase diagram agrees with our previous results: 
At zero chemical potential we have the transition point 
Wh = 0.2516. At zero temperature we have the transition 
point at Jl = 0.708, which corresponds to d = 0.3197. We 
also identify Jl = 0.470 as the position of the second order 



transition to a meson melted phase with non-zero d and 
chiral condensate c. 

The dotted green line is the line along which d = and 
corresponds to the second order meson melting transition 
from a Minkowski embedding to a black hole embedding. 
The transition generates density continuously from zero. 
The quark condensate also smoothly decreases from its 
constant value on the Minkowski embedding. We dis- 
play the continuous behaviour of the quark condensate 
across the transition in Fig 6. Note this means that the 
slope of the embedding at the UV boundary is continu- 
ous through the transition even though the embedding 
in the IR is discontinuous and topology changing. Again 
we have checked the smoothness of these parameters nu- 
merically to better than the 1% level. 

The blue line corresponds to a first order transition 
and the red dotted line is a second order transition in 
density, chiral condensate etc. The red dotted line is 
rather special in that this is a phase boundary only at 
fh = 0. This is because this phase boundary is related 
to the spontaneous breaking of chiral symmetry which 
only exists at m = 0. At finite m it must be a cross over 
region as we will discuss further in section VII 

The diagram then displays two tri-critical points. It is 
straightforward to identify where the points lie numer- 
ically. The chiral symmetry tri-critical point where the 
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FIG. 5: The phase diagram of the Ji = 2 gauge theory with a 
magnetic field. The temperature is controlled by the param- 
eter wh and chemical potential by Jl. (Parameters are scaled 
or B — 1/2R 2 in terms of parameters without tilde.) 



FIG. 7: The phase diagram of the N = 2 gauge theory with 
a magnetic field. The temperature is controlled by the pa- 
rameter wh and the density by d. (Parameters are scaled or 
B — 1/2R 2 in terms of parameters without tilde.) 
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first and second order chiral symmetry restoration transi- 
tions join lies at the point (J1,wh) = (0.267,0.201). The 
second tri-critial point where the meson melting transi- 
tions join is at (J1,wh) = (0.129,0.236). 



VI. THE PHASE DIAGRAM IN THE 
CANONICAL ENSEMBLE 

We can study the phase diagram also in the canonical 
ensemble. It is shown in Fig 7 and has the same informa- 
tion as Fig 5. The pale green region in Fig 5 lies in the 
green dotted line along the wh axis of Fig 7. The chiral 
symmetry breaking region enclosed by the red, green and 
blue lines in each figure map onto each other. Similarly 
the high temperature and density region to the upper 
right of all the lines in both plots map onto each other. 
The two double blue lines and the area between them in 
Fig 7 correspond to the single blue line in Fig 5, which 
is natural since the blue line in Fig 5 is a first order 
transition line and the density change is discontinuous. 
Thus the gray region in Fig 7 is an unstable density re- 
gion which hides in the phase boundary in Fig 5. That 
region may only be reached by super-cooling or super- 
heating since it is unstable. The true ground state at 
those densities and temperatures should be a mixture of 
the black hole and Minkowski embedding in analogy with 
the liquid-gas mixture between the phase transition's of 
water [16]. It's not clear how to realize that mixture in 
a holographic set-up. 



VII. FINITE MASS 

We next describe the evolution of the phase diagram 
with quark mass. If we move away from zero quark mass 
then the second order chiral symmetry restoration phase 
transition at T=0 but growing chemical potential be- 
comes a cross over transition. This can be seen in Fig 
2 where for ra/0 the non-zero value of the condensate 
can be seen to change smoothly with changing Jl and 
there is no jump in any order parameter. The (chiral) 
tri-critical point becomes a critical point. However, the 
other transition lines survive the introduction of a quark 
mass. 

In Fig 8. we plot the phase diagram for various quark 
mass, to, at constant B. The colors represent different 
quark masses - to — 0, 1, 1.5, 2, 3 from bottom to top are 
black, red, orange, green, and blue. The solid lines are for 
finite, fixed B. To show the influence of the magnetic field 
we also display the B = solution as the dotted lines. 
The gray line shows the motion of the critical points. 

In general the magnetic field shifts the transition line 
up and right, meaning that the magnetic field makes the 
meson more stable against the temperature/density me- 
son dissociation effect. This is important at small fh but 
negligible at large fh as expected. 



"7/ 




(a) The curves correspond to 
fh = O(black), l(red), 1.5(orange), 2(green), 3(blue) from 
bottom to top. The gray line is the path of the critical 
points. 

0.760. 



0.759 




(b)A zoom into figure (a) to show the critical point structure, 
at small chemical potential. Here we show the detail of 
in = 1 case, but a similar structure exists for every case. 

FIG. 8: The phase diagram at finite current quark mass with 
finite B (solid lines) and zero B (dotted lines). 



Both critical points survive the introduction of a finite 
to, even though it looks like there is no critical point in 
Fig 8a. Zooming in on the appropriate region at small 
chemical potential reveals the two critical point structure 
as shown in Fig 8b. Their positions, as fh changes, are 
marked by the gray line in Fig 8a. The one line represents 
the two critical points which are indistinguishable close 
on the scale of Fig 8a. The chiral symmetry critical point 
moves very close to the other critical point even for a very 
small mass (fh ~ 0.01). The interpretation of the critical 
points and the phase boundaries are the same as in the 
???. = case in the previous section. 

Notice that the black hole to black hole transition ex- 
ists even in the B=0 case as shown in Fig 8b (Right), so it 
is not purely due to the magnetic field. Nevertheless this 
transition seems not to have been reported in the previ- 
ous works [14, 16]. We believe that this is because the 
transition line between the two critical points is too small 
to be resolved on the scale of Fig 8a, which agrees qual- 
itatively with the figures in [16]. In order to find those 
transitions we had to slice the temperature down to or- 
der 10~ 3 as shown on the vertical axis in Fig 8b(Right). 
Any coarser graining would miss it. 

The final surprise relative to the previous work is that 
the meson melting transition below the critical point ap- 
pears second order in our work even in the infinite mass 
limit. To emphasize this we show a number of plots in 
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0.4. At /i ~ 0.4467 there is a Minkowski to black hole embedding transition, which is second order in both chiral 

condensation and density. 
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(b) wh = 0.7575. There are two transitions. The first (fi ~ 0.0341) is a Minkowski to black hole transition and second order in 
condensation and density. The second (p, ~ 0.03455) is a black hole to black hole transition and first order. 
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0.7587. At fi ~ 0.0321 there is a Minkowski to black hole embedding transition, which is first order in both chiral 

condensation and density. 
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(d) w H = 0.762. At p > 

wh = 0.7658 
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' 0.0235 there is a Minkowski to black hole embedding transition, which is first order in both chiral 
condensation and density. 





(c) Above wh = 0.7658 only a black hole embedding (Red) is stable configuration. 



FIG. 9: Chiral condensation (Left), density (Middle), and the grand potentials (Right) for massive quarks (m = 1) at B 
at a variety of temperatures that represent slices through the phase diagram Fig 8a. 
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the B = theory in Fig 9. 

Since the scaled variables (16) cannot be used at B = 
0, (20) and (18) read in terms of the original coordinates: 



fl(wH,fJ,) : 



-5 



N f T D7 Vo\ 

4 4 

dp .... H 



(1 + (L') 2 ) (w* + w%\ „6 

p •> 

(21) 



where 



/' 
K 



dp d 



PH ^ + W % 

4 i 4 \ 2 
W + W 



w*-w% / ! + (£') 



*1 P e 



K 

,4 



(w 4 + Wjj) 



(22) 
(23) 



p := \pi ira'A t (oo) , d 



7 d (24) 



By the same procedures as in the previous sections we 
get Fig 9. Compared to Fig 3, the left column of Fig 9 
is the chiral condensate instead of the embedding con- 
figurations. In Fig 3 there is always a red black hole 
embedding, which corresponds to the flat embedding at 
zero quark mass. It is not present at finite quark mass. 

At very low temperature the transition is Minkowski to 
black hole and second order in the condensate and den- 
sity (Fig 9a) . As the temperature goes up a new black hole 
to black hole transition pops up by developing a 'swal- 
low tail' in the grand potential - this transition is first 
order in the condensate and density (Fig 9b). As temper- 
ature rises the 'swallow tail' grows continually and even- 
tually "swallows" the second order Minkowski to black 
hole transition (Fig 9 c,d). i.e. At higher temperature the 
second order Minkowski to black hole transition enters 
an unstable regime and plays no role any more. Instead 
only the first order Minkowski to black hole transition is 
manifest. Finally the Minkowski embedding becomes un- 
stable compared to the black hole embedding(Fig 9c). At 
an even higher temperature the Minkowski embedding is 
not allowed and only a black hole embedding is available 
(Not shown in Fig 9). 

These results all match with our work at finite B and 
increasing mass, confirming those results and our phase 
diagrams already presented. 



VIII. COMPARISON TO QCD 

We have computed the phase diagram for a particular 
gauge theory using holographic techniques. There are 
many differences between our theory and QCD: the the- 
ory has super partners of the quarks and glue present; it 
is at large N and small Nf, so quenched (and we have 
only computed for degenerate quarks to avoid complica- 
tions involving the non-abelian DBI action); the theory 























\ 









(a)Standard scenario 



(b)Exotic Scenario 



FIG. 10: Two possible phase diagrams for QCD with the 
observed quark masses, (a) is the standard scenario found 
in most of the literature but a diagram as different as (b) 
remains potentially possible according to the work in [3]. We 
have not included any color superconducting phase here at 
large chemical potential. 



has deconfincd glue for all non-zero temperature; the the- 
ory has a distinct meson melting transition. In spite of 
these differences the phase diagram for the chiral con- 
densate shows many of the aspects of the QCD phase 
diagram so we will briefly make a comparison here. 

The QCD phase diagram is in fact not perfectly 
mapped out since there have only recently been lattice 
computations attempting to address finite density [3]. 
The phase structure also depends on the relative masses 
of the up, down and strange quarks. The standard the- 
oretical picture [1-3] for physical QCD is shown in Fig 
10a. At zero chemical potential the transition with tem- 
perature is second order (or a cross over with massive 
quarks). At zero temperature there is a first order tran- 
sition with increasing chemical potential (ignoring any 
superconducting phase) . These transitions are joined by 
a critical point. Comparing to our theory in Fig 5 we see 
that the transitions' orders are reversed and the pictures 
look rather different. 

In fact though as argued in [3] the picture could be very 
different in QCD. At zero quark mass the finite temper- 
ature transition is first order and whether it has changed 
to second order depends crucially on the precise physical 
quark masses. Similarly whether the finite density tran- 
sition is truly first order or second order depends on the 
exact physical point in the m Ut d, m s , P, T volume. Ar- 
guments can even be made for a phase diagram matching 
that in Fig 10b which then matches the structure of the 
chiral symmetry restoring phase diagram of the theory 
we have studied. For the true answer in QCD we must 
wait on lattice developments. Clearly our model will not 
match QCD's phase diagram point by point in m Ul d, m s , 
fi, T volume but it provides an environment in which 
clear computation is possible for structures that match 
some points in that phase space. 

Finally, we note a more general point that seems to 
emerge from the analysis. The introduction of a chemi- 
cal potential weakens the first order nature of the tran- 
sitions in our analysis. This matches with results found 
in QCD on the lattice. The weakening of the first order 
phase transition is demonstrated for the chiral transition 
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in the light quark mass regime [35, 36], and is shown for 
the deconfincnient transition in the heavy quark mass 
regime [37, 38]. 
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